Fishtail effect and vortex dynamics in LiFeAs single crystals 
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We investigate the fishtail effect, critical current density ( J c ) and vortex dynamics in LiFeAs single 
crystals. The sample exhibits a second peak (SP) in the magnetization loop only with the field || 
c-axis. We calculate a reasonably high J c , however, values are lower than in 'Ba-122' and 'llll'-type 
FeAs-compounds. Magnetic relaxation data imply a strong pinning which appears not to be due to 
conventional defects. Instead, its behavior is similar to that of the triplet superconductor Sr2Ru04. 
Our data suggest that the origin of the SP may be related to a vortex lattice phase transition. We 
have constructed the vortex phase diagram for LiFeAs on the field-temperature plane. 

PACS numbers: 74.25. Ha, 74.25. Sv, 74.25.Wx 



The fishtail effect or the anomalous second peak (SP) 
in field (if) dependent magnetization (M) loops has been 
a subject of intense research topic in the field of super- 
conductors with both low and high transition temper- 
atures (fc)»i~— This phenomenon is realized with the 
enhanced irreversibility in isothermal M(H) or, equiva- 
lently, the enhanced critical current density ( J c ) at high 
fields apart from the central peak which occurs around 
zero field. In type-II superconductors, the magnetic fields 
above the lower critical field (if c i) penetrate the bulk of 
the sample in the form of flux lines or vortices. Differ- 
ent mechanisms based on the vortex dynamics have been 
discussed to explain the SP in high-Tc cuprates which in- 
clude inhomogeneity of the sample^ matching effect, 4 
surface barriers^ geometrical effects^ dynamic effects, 7 
structural phase transition in the vortex lattice (VL)^ 
vortex order-disorder phase transition,— ~— crossover from 
elastic to plastic creep^ etc. However, in spite of plenty 
of studies dealing with this phenomenon, the general un- 
derstanding lacks a converging trend and the proposed 
models appear to be more sample specific. 

The recent discovery of superconductivity (SC) in Fe- 
based pnictides^ has renewed the interest in vortex dy- 
namics. Similar to cuprates, pnictides are also layer- 
based superconductors, and exhibit a high Tq and type-II 
nature. In contrast, pnictides have less anisotropy and 
larger coherence length (£), thus raising question how 
these influence the vortex dynamics in these materials. 
The appearance of a SP in M(H) is not an universal phe- 
nomenon in different families of pnictides. For the '122' 
family (AFe2As2, A = Ba, Sr, Ca, etc.) the appearance 
of a SP is sensitive enough to the compositional elements. 
The pronounced SP has been observed, for example, in 
both hole- and electron-doped Ba-122 compounds, which 
has been ascribed to various mechanisms but it re- 
mains absent in doped Ca-122 compounds. 22 Similarly, a 
SP is not consistently seen in the '1111' (REFeAsO, RE = 
La, Nd, Ce, Sm, etc.) and in the '11' (FeTe) families 
However, to the best of our knowledge the vortex dynam- 
ics have not been studied in the '111' family (AFeAs, A 
— Li, Na, etc) up to now. 

Here we study the fishtail effect and the vortex dy- 
namics in single crystals of LiFeAs which belong to the 



111 family of pnictides. LiFeAs is an oxygen free com- 
pound where superconducting active FeAs layers are 
separated by Li atoms along the c-axis^ Remarkably, 
LiFeAs exhibits SC in absence of any notable Fermi sur- 
face nesting and static magnetism^ 7 - however, the pres- 
ence of antiferromagnetic fluctuations in the normal state 
is inferred from the nuclear magnetic resonance (NMR) 
measurements.— Interestingly, recent experimental NMR 
results 2 ^ and also theoretical calculation s 30 ' 31 indicate a 
possible p-wave SC state in LiFeAs, which is significantly 
different from other families within pnictides. As LiFeAs 
is nonmagnetic and does not require any chemical doping 
to become superconducting, therefore the FeAs layers are 
more homogeneous and the crystal is devoid of coexist- 
ing magnetic phases which make it an ideal playground 
to study the vortex dynamics. Based on the sharpness 
of the rocking curve, recent small angle neutron scatter- 
ing (SANS) measurements have revealed that that VL 
in LiFeAs exhibits no long range order, however, better 
ordering than the doped Ba-122 compounds i 18 ' 20 ' 32 

We have investigated the properties of the vortex state 
in LiFeAs by means of isothermal M(if) and magnetic 
relaxation measurements which are the most extensively 
used tools for a variety of superconducting materials! 33 ' 34 
Our results imply a pronounced SP in both M(if) and 
J c (ff) at low temperatures (T) with the applied field par- 
allel to the c-axis. We do not find a SP with if | \ab plane. 
We determine the J c which is reasonably high, however, 
lower than those for doped Ba-122 and 1111 compounds. 
The magnetic relaxation data, on the other hand, imply 
a nonlogarithmic time (t) dependence with the estimated 
normalized magnetic relaxation rate (S) being very low 
(even lower than in Ba-122 and 1111 compounds) signi- 
fying a strong pinning for LiFeAs. This specific behavior 
of the magnetic relaxation is similar to those of triplet 
superconductors, i.e., S^RuC^i 3 -^ Our data also indicate 
that the SP in LiFeAs may be due to a VL phase tran- 
sition. From the M(H) plots we have constructed the 
vortex phase diagram on the ff-T plane for LiFeAs. 

Single crystals of LiFeAs have been grown using the 
self-flux method as detailed in Ref. [36J The good 
quality and homogeneity of the crystals are confirmed 
by a exceptionally small nuclear quadrupole resonance 
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FIG. 1: Temperature dependence of Xvoi for LiFeAs deduced 
from dc-M measured following ZFC and FC protocols. The 
data have been collected in H\\c — 20 Oe. The inset shows M 
vs H plots measured at 10 K with the field parallel to both 
c-axis and afc-plane. 



(NQR) linewidth of 64 kHz, a very low residual resis- 
tance of 0.025 mficm and a sharp transition in specific 
heat data i 36 ' 37 For the present studies, two crystals (5*1 
and S2) of the same batch with rectangular shape have 
been selected. For the magnetic hysteresis loop sample 
SI (3.53 x 2.5 x 0.21 mm 3 ) and for the magnetic relax- 
ation measurements sample 52 (2.89 x 2.16 x 0.38 mm 3 ) 
have been used. Magnetization have been measured in a 
Quantum Design MPMS-XL SQUID. Adequate care has 
been taken to avoid the exposure of the sample to air be- 
fore mounting it in the magnetometer. All the M(H) and 
M (t) measurements have been performed after cooling 
the sample in zero magnetic field from much above Tq to 
the specific temperature. The M{H) loops have been in- 
vestigated with the field up to 50 kOe. For the relaxation 
measurements, the magnetization has been measured as 
function of time for about 8000 s. 

The main panel of Fig. 1 presents the temperature de- 
pendence of the volume susceptibility (xvoi) measured 
following the zero-field-cooled (ZFC) and field-cooled 
(FC) protocols for magnetization measurements. Xvoi 
has been deduced from the measured dc-M with a field of 
20 Oe applied parallel to the crystallographic c-axis. The 
data has been corrected for demagnetization effects £2- It 
is evident from Fig. 1 that the sample exhibits bulk SC 
as characterized by the diamagnetic signal at low T. The 
sharp transition as well as the high value of Xvoi in Mz fc 
demonstrate the high quality of our crystal. We deter- 
mine Tq from the bifurcation point between ZFC and 
FC branches of the magnetization to be around 16.5(5) 
K. In the inset of the Fig. 1 we have plotted the M(H) 
data at 10 K with H parallel to both c-axis and a&-plane. 
With increasing H for H\\c, the magnetic irreversibility 
(Mi rr ) initially decreases showing a minimum at a field 
H m . On further increase in H, Mi rr increases and ex- 
hibits a peak (SP) at a field H p . However, we do not find 
any trace of a SP for H\ |a6-plane. While this significant 
anisotropic behavior in appearance of the SP is similar to 
other pnictide superconductors^ but it remains different 
from cuprates, i.e., Lai_ x Sr ;E Cu04ii2 

Since the SP in M(H) is only evident for H\\c in this 
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FIG. 2: (Color online) (a) The isothermal M vs H loops 
recorded at different temperatures with H\\c for LiFeAs. (b) 
The critical current J c as calculated from the M(H) loops in 
(a) as function of field for different temperatures; for details 
see text. 



compound, we have collected the M (H\\C) isotherms at 
different T within the SC regime in order to understand 
the SP characteristics. The data are plotted in Fig. 2a 
where M vs H loops are quite symmetric with respect to 
both the sweeping direction as well as the polarity of the 
magnetic field. At low T, however, even though the onset 
of the SP is evident, the SP cannot be observed within 
the measurable field range. Interestingly, M(H) at 2 K 
exhibits irregular jumps close to H = 0. These jumps 
are commonly known as 'flux jump' effects^ and will be 
presented elsewhere in more detail. With the increase 
in T, a clear SP can be observed in the M(H) loops. 
Moreover, we find a field Hi rr above which Mi rr in data 
vanishes. With increasing T, all the characteristics fields, 
i.e., H mi H p and Hi rr decrease, and their T variation will 
be discussed in a later section. 

From the magnetic irreversibility in M(H) we have 
calculated the critical current J c exploiting the Bean's 
critical state model^ J c — 20AAf/[a (1 — a/36)], where 
AM = Mdn - M up , M up and Mdn are the magnetization 
measured with increasing and decreasing field, respec- 
tively, and a and b (b > a) are the dimensions of the 
crystal perpendicular to the applied H . The unit of AM 
is in emu/cm , a and b are in cm and the calculated J c is 
in A/cm 2 . The calculated J C (H) has been plotted in Fig. 
2b for different T. The variation in J C (H) is nonmono- 
tonic and exhibits a broad peak (SP) in the high field 
region, which is in line with AM in Fig. 2a. At low T, 
J c is rather high, however, its value still being lower than 
in doped Ba-122 and 1111 compounds where J c ~ 10 6 or 
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FIG. 3: (Color online) (a) The energy barrier Uo/ks obtained 
by fitting of Eq. f as a function of temperature, (b) The 
magnetic relaxation rate S as function of temperature in H\\c 
= 10 kOe as calculated using Eq. 2 at t = 1000 s. The inset 
exemplary shows the best fit of the magnetic relaxation data 
at T = 10 K and H\\c = 10 kOe using Eq. 1. 



even higher i^ r 17 i 19 ' 24 This is in agreement with the level 
of disorder as revealed from the SANS studiesj 18 ' 20 ' 32 

To understand the origin of the SP and the reasonably 
high J c - values in LiFeAs, we have studied detailed vor- 
tex dynamics in this compound by means of T- and in- 
dependent magnetic relaxation measurements! 33 ' 34 Mag- 
netic relaxation in superconductors is a result of non- 
equilibrium spatial arrangement of vortices due to pin- 
ning sites. External applied magnetic field exerts Lorentz 
forces on the vortices resulting in their movement, which 
causes a change in M(t). In contrast to the original 
Anderson-Kim model jsi our relaxation data exhibit a 
nonlogarithmic time dependence and can be best fitted 
with the interpolation formula: 33 



M(t) = M a 



1 



nk B T 



In 



-l/p 



(1) 



where ks is the Boltzmann constant, Uq is the energy 
barrier height in absence of a driving force, to is the char- 
acteristic relaxation time (usually ~ 10 _6 s for type-II su- 
perconductors), and /j, is the field-temperature dependent 
parameter. This formula yields the normalized magnetic 
relaxation rate S[= (l/M)dM/d\n(t)] as^ 



S(t) 



fc R T 



U Q + [ik B Tln(t/t ) 



(2) 



Magnetic relaxation has been measured at different T 
in H\\c = 10 kOe. We find a very slow relaxation, i.e., 
there is only a 4% change in magnetic moment at H = 
10 kOe and T = 10 K within a time period of 8000 s 
which is much lower than that observed in cuprates and 
122-pnictides^i^ 3 . U /k B (T) and S(T) as extracted from 
the fitting of the data exploiting Eqs. 1 and 2 are plot- 
ted in Fig. 3a and 3b, respectively. One representative 
fitting of our data has been included as an inset in Fig. 
3b. Surprisingly, we find a very high value of to of the 
order of 10 s, which is orders of magnitude higher than 
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FIG. 4: (Color online) (a) The field dependence of Uo/ks 
and S extracted from the magnetic relaxation data at 10 K 
exploiting Eqs. 1 and 2. The vertical dashed lines mark H m 
and H p ; details see text, (b) The same physical properties 
Uo/ks and S shown for 5 K. 



for other families in pnictides and cuprates i 15 ' 22 ' 33 ' 34 For 
LiFeAs the value for the energy barrier Uo / k B increases 
with T, however, above 10 K it drops down due to larger 
thermal fluctuations. A similar behavior can be detected 
for S(T), which also shows a non- monotonous behavior 
(calculated at t = 1000 s) with a steep decrease followed 
by an increase upon lowering temperature. 

In addition to the temperature dependent relaxation 
studies, M(t) has been measured in different fields (if||c) 
along the hysteresis loop at constant T. Similarly, we 
have extracted Uo/ks and S using Eqs. 1 and 2. Fig. 
4a presents the results for Uo/ksiH) and S(H) at 10 K. 
Here, the vertical dashed lines in Fig. 4a represent the 
minimum and maximum (SP) in M(H) at 10 K (sec Fig. 
1 and 2). It is evident from the plot that Uo/ks initially 
increases with applied fields and then decreases, show- 
ing a peak at a field in between H m and H p , where the 
related calculated S(H) exhibits a minimum. This varia- 
tion in S(H) is similar to doped Ba-122 compounds . 15 ' 21 
To examine the behavior of Uo/ksiH) and S(H) at low- 
T where the SP shifts significantly to higher field values, 
we have calculated the parameters at 5 K following the 
same method. Our results show that although the field, 
where Uo/ks {S) exhibits a maximum (minimum), in- 
creases from wlO kOe at 10 K to «20 kOe at 5 K (Fig. 
4b), this variation does not scale with the large increase 
in Hp at low T (see Fig. 2a). Note, that the calculated 
S in Figs. 3 and 4 is very low, and that the values even 
go below those for Ba-122 and 1111 compounds, where 
S > O.Ol^i^i^ 4 The low value of S paired with our 
very high io-value imply a high pinning in LiFeAs. This 
is quite intriguing as J c in LiFeAs is lower than in Ba- 
122 and 1111 materials; 1 ^ 17 ' 19 ' 24 which means that other 
type of pinnings rather than the conventional defects are 
active in LiFeAs. Indeed, our crystal is of good qual- 
ity and homogeneous as discussed earlier. However, such 
slow magnetic relaxation has been observed in the p-type 
superconductor Sr 2 Ru04 (S ~ 10~ 3 ), where it has been 
pointed out that the observed strong pinning is related to 
the superconducting phase of triplet nature rather than 
conventional defects^ This raises the question if the ob- 
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FIG. 5: (Color online) The vortex phase diagram on the H-T 
plane with H\\c for LiFeAs. Solid lines represent fits of the 
data using the functional form H X (T) = H X (0)(1 - T/T c ) n 
(see text). 

served low S- value in combination with the observed J c 
in our single crystal of LiFeAs might be within the lines 
of other studies, which recently proposed triplet pairing 
m LiFeAs^S^i but this requires further detailed investi- 
gations in the future. 

From the so far obtained data we have constructed the 
vortex phase diagram on the H-T plane for LiFeAs (Fig. 
5). Above Hi rr , vortices are in a liquid state. Below Hi rr , 
they are in a solid state, however, its nature changes 
between different field regimes indicated in the figure 
(named I, II and III). All the characteristic fields show 
a strong T dependence where the data can be fitted well 
with the functional form H x (T) = H x (0)(l-T/T c ) n . We 
obtain H irr (0) = 299.5(9) kOe and n = 1.46(4), H p (0) = 
105.7(8) kOe and n = 1.66(3), H m (0) = 29.1(4) kOe and 
n = 1.56(3). The values of the exponent n are reasonably 
consistent with those for other pnictidesi 15 i 16 ' 23 ' 24 

Now we discuss the origin of the SP effect and the 
intriguing vortex dynamics in LiFeAs. The strong T de- 
pendence of the transition lines (Fig. 5) discards the 
possibility that the SP in LiFeAs arises due to a vor- 
tex order-disorder phase transition^ - — Moreover, an ab- 
sence of 'mirror-image' correlation between M(H) and 
S(H) in Figs. 2 and 4 imply that the dynamic model 7 
is not valid in the present case. Although, the nature 
of Uo/kB(H) in Fig. 4 has qualitative similarity to the 
model which predicts the SP being associated with a 
crossover in flux dynamics from elastic (< H p ) to plas- 
tic (> Hp) creep with increasing field^ the peak in 



Uo/ks(H) for LiFeAs occurs much below the SP, and 
in low-T at 5 K this mismatch increases. Moreover, this 
model predicts H p oc [1 — (T/Tc) 4 ] 1 ' 4 which does not de- 
scribe the present H p (T) dependence shown in Fig. 5. 
Henceforth, the validity of this model for LiFeAs is ques- 
tionable. 

On the other hand, a structural phase transition in 
the VL is also an attractive and possible model which 
argues that the SP is associated with the transforma- 
tion of a hexagonal VL to a square one with field.— The 
square structure is supported by the fourfold symmetry 
of the intervortex interaction which can originate in var- 
ious situations, like, for the anisotropic (d-wave) nature 
of the superconducting gap as in La2_ x Sr. r Cu04( 4 ^ for 
materials with low Ginzburg-Landau (GL) parameters 
(k) as in YM2B2C,— in the extended GL theories with 
more than one order parameter as for the p-wave super- 
conductor Sr2Ru04, 44 etc. Recently, such scenario of a 
structural phase transition in the VL has been proposed 
in Ba(Feo.925Coo.o75)2As2 where a minimum in S(T) and 
S(H) has been found) 21 - Altogether, taking into account 
a similar behavior with a minimum in S(T) and S(H) 
in LiFeAs, paired with a comparatively low value of k 
(«30)— and the proposed p-wave SG 2 ^— in LiFeAs, a 
structural phase transition in the VL is a possible sce- 
nario for the existence of the SP in this compound. How- 
ever, further investigations including microscopic probes 
are required to confirm these observations. 

In summary, single crystalline LiFeAs exhibits a SP in 
the M(H) loop with iJ||c-axis. The calculated J c s are 
reasonably high, however, the values are lower than in 
the doped Ba-122 and 1111 compounds. We find an ex- 
traordinary slow magnetic relaxation implying a strong 
pinning which appears not to be related to conventional 
defects. Instead, the behavior of magnetic relaxation is 
similar to the p-wave superconductor S^RuO-j. We have 
constructed the vortex phase diagram on the H-T plane 
for LiFeAs, with the characterized fields Hi rr , H p and 
H m showing a strong T dependence. In accordance with 
recent investigations on Ba(Feo.925Coo.o75)2As2 our data 
imply that the SP in LiFeAs most likely originates from 
a VL phase transition. Nonetheless, further studies in- 
volving microscopic probes are required to comprehend 
the SP and vortex dynamics in this compound. 
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